Improved constraints on non-Newtonian forces at 10 microns 
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Several recent theories suggest that light moduli or particles in "large" extra dimensions could me- 
diate macroscopic forces exceeding gravitational strength at length scales below a millimeter. Such 
new forces can be parameterized as a Yukawa-type correction to the Newtonian potential of strength 
q relative to gravity and range A. To extend the search for such new physics we have improved 
our apparatus utilizing cryogenic micro-cantilevers capable of measuring attonewton forces, which 
now includes a switchable magnetic force for calibration. Our most recent experimental constraints 
on Yukawa-type deviations from Newtonian gravity are more than three times as stringent as our 
previously published results, and represent the best bound in the range of 5 — 15 /im, with a 95% 
confidence exclusion of forces with \a\ > 14,000 at A = 10 fim. 

PACS numbers: 04.80.Cc 



I. INTRODUCTION 

A number of theories of physics beyond the Standard 
Model suggest that new physics related to gravity may 
appear at sub-millimeter length scales. For example, 
light moduli from string theory ELB or exotic particles 
in "large" extra dimensions 0, [j, Q can mediate forces 
exceeding gravitational strength that can be observed in 
a tabletop experiment. Such phenomena can be param- 
eterized as a Yukawa-type correction to the Newtonian 
potential of strength relative to gravity a and range A. 
For two masses mi and rri2 separated by distance r, the 
gravitational potential is modified to 
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To search for such new forces, we have improved our cryo- 
genic apparatus utilizing silicon micro-cantilevers with 
attonewton force sensitivity [f| 0], which now includes 
a magnetic method for force calibration. The cantilever 
is loaded with a rectangular gold prism fabricated by 
focused-ion-beam milling that serves as a test mass for 
the experiment. The driving (source) mass is moved 
horizontally beneath the cantilever at a nominal vertical 
face-to-face separation of 25 fim. The force between the 
masses is deduced from the displacement of the cantilever 
as measured by a fiber-coupled laser interferometer. We 
perform the measurement at the cantilever resonance fre- 
quency, typically of order 300 Hz, while the mechanical 
driving motion occurs at a sub-harmonic, typically one- 
third. This is achieved by implementing a density mod- 
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ulation in the drive mass, consisting of alternating gold 
and silicon sections. 

For the magnetic calibration, a Co/Pt multi-layer film 
is deposited on the test mass. The permanent mag- 
netic moment couples to a magnetic field gradient pro- 
duced by current flowing across the meandering gold sec- 
tions in the drive mass device. The current is turned 
off for the Yukawa-force search. A shield composed of 
high-magnetic-permeability material encloses the cryo- 
stat to prevent the Earth's field from magnetizing the 
drive mass. The amplitude and phase of a magnetic or 
Yukawa signal will change in a predictable way as we vary 
the equilibrium position of the drive mass oscillation. We 
utilize this scanning technique as an additional handle to 
distinguish a signal from background forces. In this pa- 
per, the newest data and error analysis are described, and 
the latest experimental constraints on Yukawa-type de- 
viations from Newtonian gravity are presented. Finally, 
future directions are discussed. 



II. EXPERIMENTAL SETUP 

Many of the details of the experimental probe and vac- 
uum cryostat are described in previous work 0, @, [!]• 
The single-crystal silicon cantilevers are 250 fim long, 50 
fim wide, 0.3 /xm thick, and have a spring constant k of 
approximately 0.0062 N/m. The thermal-noise-limited 
force sensitivity at 10 K is approximately 200 aN/Hz 1 / 2 
for a typical low-temperature quality factor Q of 80,000. 
The displacement of the cantilever beam is measured by 
using a fiber-coupled laser interferometer where a Fabry- 
Perot cavity is created between the end of the fiber op- 
tic and the cantilever loaded with a test mass 0. The 
force F on the cantilever is deduced from the measured 
displacement z, which is enhanced on-resonance by the 
large quality factor: z = FQ/k. The 1.5 /ig test masses 
measure 54 x 54 x 27 /im 3 and are cut from the edge of 
a 27 /zm thick gold foil using a 20 nA gallium focused- 
ion-beam (FIB). This technique produced masses with 



more regular shapes than those used in previous versions 
of the experiment [7[ • A SEM image of a FIB-cut mass 
attached to a cantilever is shown in the inset of Fig. [U 
The improvement allows a flatter surface to be presented 
to the drive mass, allowing more force sensitivity, and 
results in an improved interferometer signal. The driv- 
ing (source) mass is mounted on a piezoelectric bimorph 
which actuates in the y-direction (as shown in Fig. [T]) 
beneath the cantilever with an amplitude of ~ 120 /im 
about its equilibrium position. The drive mass consists of 
alternating 100 /im-wide bars of gold and silicon that are 
approximately 100 /im deep and 1 mm long. As the bi- 
morph oscillates with a sinusoidal motion, a time- varying 
gravitational force is exerted on the cantilever which can 
occur at the driving frequency of the bimorph as well as 
at higher harmonics, depending on the drive amplitude 
and spatial orientation of the test and drive masses. Sim- 
ulation indicates a maximal gravitational coupling at the 
third harmonic of the drive frequency for a bimorph am- 
plitude of 133 ;um, with only a slight reduction (~ 10%) 
at 120 /im 10]. A gold-coated silicon nitride shield mem- 
brane separates the cantilever from the drive-mass and 
provides attenuation of electrostatic and Casimir back- 
ground forces from the oscillating drive mass. The drive 
mass is covered with a smooth 1 /tm-thick plane of sil- 
icon, followed by aluminum oxide and gold to further 
suppress modulations in electrostatic or Casimir forces 
associated with the sections of alternating density. A 
schematic (not to scale) showing the cantilever, test and 
drive masses, shield, and piezoelectric bimorph actuator 
is shown in Fig. [TJ A piezoelectric stack actuator is also 
included near the base of the cantilever to allow tests 
with deliberate excitation and to facilitate interferome- 
ter offset control and characterization. The voltage from 
the signal photodiode in the interferometer and the volt- 
age applied to the bimorph actuator are recorded at 10 
kHz on a data acquisition (DAQ) device connected to a 
PC. The cantilever signal (at the third harmonic of the 
drive frequency) is averaged with respect to the phase of 
the drive signal. A series of time records is collected for 
each data point. An FFT of the interferometer data is 
performed to determine the amplitude and phase at the 
third harmonic of the drive signal. 

The amplitude and phase of either a magnetic or 
Yukawa signal will change in a predictable way as we 
vary the y-equilibrium position of the drive mass oscilla- 
tion in situ. The x,y, and z-separation between the drive 
mass and test mass is controlled by a partially (y,z) mo- 
torized manipulator stage at the top of the cryostat, and 
the tilt and x,y, and z-separation between the masses 
are detected by capacitive sensors @, 0] ■ The expected 
Yukawa signal for a typical vertical mass separation and 
bimorph amplitude is shown in Fig. [5J The vertical axis 
shows the calculated amplitude and phase of a Yukawa 
(a = 1000, A = 18 /im) force at the third harmonic of 
the drive (which equals the resonance frequency of the 
cantilever). The horizontal axis shows the y-equilibrium 
position of the oscillation. Note that the third harmonic 
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FIG. 1: (color online). A schematic (not to scale) showing 
the cantilever, test and drive masses, shield, and piezoelec- 
tric bimorph actuator. Figure is adapted from Ref. (Tl| . 
Not shown are the metallization of the shield and drive mass 
ground plane. A Fabry-Perot cavity is formed between the 
bottom of the optical fiber and top of the test mass and is 
used to interferometrically measure the cantilever displace- 
ment. The coordinate axes shown in the lower inset are used 
throughout the paper to describe the orientation of the drive 
mass and test mass. Upper inset: SEM micrograph of FIB- 
fabricated test mass cut from 27 /tm thick gold foil attached 
to cantilever. 



of the force (along with other odd harmonics) vanishes 
when the equilibrium position of the drive mass is such 
that a gold bar or silicon bar is centered underneath the 
cantilever. Conversely, the force is maximized when the 
interface of the bars is centered beneath the cantilever. 
This results in a 100-/im spatial periodicity for the am- 
plitude of a Yukawa-type signal. 



III. MAGNETIC CALIBRATION 

Having a magnetic force for calibration that is switch- 
able in situ allows a great improvement in systematics 
by providing additional information about the expected 
phase structure of a Yukawa-type signal and by allow- 
ing multiple days of data to be analyzed together, as 
explained in Sec. |V1 When current is applied in the 
meandering pattern of gold bars in the drive mass, the 
resulting magnetic held gradient above the pattern cou- 
ples to a magnetic moment on the cantilever to provide 
a measurable force. This allows the relative position of 
the test mass with respect to the drive mass to be deter- 
mined, as explained in the following subsection. When 
the current is turned off, for the Yukawa-measurement 
to have a clean background, it is necessary to eliminate 
or suppress the force due to the differing magnetic sus- 
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FIG. 2: (color online). Calculated vertical component of the 
Yukawa (a = 1000, A = 18/im) force amplitude and phase at 
the third harmonic of the drive frequency (3a>), for 119 /jm 
of bimorph amplitude and 31.5 fim edge-to-edge z-separation 
between masses with 6 /im tilt across the drive mass in the 
scanning (y) direction. The horizontal axis shows the y- 
equilibrium position of the drive mass oscillation. For ref- 
erence, the 100 /im-wide gold bars are centered in the cal- 
culation at y — —250 /im, y — —50 /im, y — 150 /mi, and 
y — 350 /im, with silicon bars in between. 



ceptibility of the gold and silicon bars in the drive mass 
when the drive mass becomes magnetized in an ambient 
magnetic field such as the Earth's (see Sec. IIIB). 

A. Principle of the calibration 

The vertical force on the cantilever with permanent 
magnetic moment m c in a magnetic field B is given by 

F z = {m c -V)B Z . (2) 

Also a torque on the cantilever f = m c x B can result in 
an effective force 

K = Tz/lc (3) 

where l c is the length of the cantilever. Simulations are 
performed to model the magnetic field from current flow- 
ing in the drive mass pattern, and the magnitude of the 
magnetic moment can be deduced from Eqs. () 12 1 - i [) and 
the measured force on the cantilever, provided the direc- 
tion of the magnetic moment is known. If the direction of 
the net moment is uncertain, only an approximate mag- 
nitude can be determined (see Fig. [3]). 

As we scan the y-cquilibrium position of the oscillation 
in situ, the magnetic signal at the third harmonic of the 
motion frequency varies with a periodicity of 200 /im, 
corresponding to the spacing between the gold bars. In 
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FIG. 3: (color online). Calculated magnetic force at the third 
harmonic of the drive frequency, per unit magnetic moment 
per ampere of current through the meander for varying y- 
equilibrium position of the drive mass oscillation. Results are 
shown for 119 fim of bimorph amplitude and 30 /im face-to- 
face z-separation, for net magnetic moment in the z,y, and x 
directions, respectively. 



particular, as the y-equilibrium position of the drive mass 
oscillation is scanned, the amplitude of the force on the 
cantilever goes through a series of maxima and minima 
with a separation of 200 /im between adjacent maxima 
and minima. Correspondingly, the phase of the force on 
the cantilever experiences a shift of tt across each mag- 
netic minimum. The particular locations of the maxima 
and minima with respect to the position of the gold bars 
depends on the direction of the net magnetic moment. 
This is in contrast to the Yukawa signal, for which the 
maxima occur when the equilibrium position of the drive 
mass oscillation is such that the test mass is centered 
over the interface between the gold and silicon bars. The 
calculated magnetic force on the cantilever per ampere 
of current in the drive mass per unit magnetic moment is 
shown in Fig. [3] for the cases of a net moment aligned in 
the x,y, or z directions, respectively. The vertical separa- 
tion between the drive and test masses is taken to be 40 
/tm and the bimorph amplitude is 119 /im. The horizon- 
tal axis shows the y-equilibrium position of the bimorph 
oscillation. 

To obtain a proof-of-principle measurement of the 
magnetic response, an early test was performed with a 
magnetic test mass consisting of a gold rectangular prism 
coated with 100 nm of nickel [7]. The y-equilibrium po- 
sition of the drive mass oscillation was varied across the 
entire drive mass pattern as the force on the cantilever 
was studied with and without current running through 
the drive mass meander. With current in the drive mass 
meander, the force occurred with a 200 /im periodicity as 
expected, corresponding to the spacing between the gold 
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bars carrying current. With no current flowing across 
the drive mass, the force occurred with a 100 fim period- 
icity and is due to the magnetization of the drive mass 
in the ambient earth's magnetic field. It is necessary to 
eliminate or suppress this susceptibility-dependent force 
in order to detect a new non-Newtonian signal. 

The phase of the magnetic force both with and without 
current flowing across the drive mass is identical modulo 
7T, as seen in earlier work Q and verified in simulation. 
Simulation indicates that the phase of any Yukawa signal 
will also be identical to the phase of the magnetic force 
modulo 7r. This is a useful handle in the identification of 
a Yukawa force (see Sec. PVT) . 



B. Magnetic shield 

In order to minimize the effect of the Earth's field, we 
employ a dual-layer magnetic shield assembly consisting 
of 0.050" thick high magnetic permeability material pro- 
vided by Amuneal Corporation. The shield is designed so 
it can be placed in position surrounding the cryostat. For 
data collection the cryostat is hung from the ceiling and 
hangs freely within the inner diameter of the magnetic 
shield. Numerical simulations performed by Amuneal 
Corporation [l2j indicated an expected transverse shield- 
ing factor of 260 and longitudinal shielding factor of 160. 
At Stanford, the earth's magnetic field is primarily ver- 
tical with a z-component of 45 ^T, so the longitudinal 
figure is relevant. The shield is degaussed prior to the 
experimental runs. 



C. Co/Pt multilayer films 

The ideal magnetic moment would have a known mag- 
nitude and direction. Cobalt/platinum multi-layer films 
have been shown to have out-of-plane anisotropy for lim- 
ited ranges of cobalt and platinum thickness [l3|, [l4[ • In 
order to obtain out-of-plane anisotropy the cobalt film 
thickness must be between 2.5 A and 11.5 A. The plat- 
inum spacing layers must also be larger than ~ 8 A [lT 
For the test masses, 3 A Co/ 15 A Pt/ 3 A Co/ 30 
Pt is sputtered on a 3 mm square gold foil substrate, on 
top of a seed layer of 100 A Pt deposited at 350 C. In- 
creasing the platinum thickness from 11 to 15 A resulted 
in improved out of plane anisotropy. The foil is then in- 
serted into an MPMS-SQUID [ljl, and a hysteresis loop 
is studied for the film magnetized both out-of-plane and 
in-plane for fields of -0.318 MA/m to 0.318 MA/m. The 
coercive field is around 0.056 MA/m and full saturation 
occurs above 0.24 MA/m. The results shown in Fig. Q] 
indicate an out-of-plane easy axis, although the in-plane 
magnetization loop also shows significant hysteresis. We 
attribute this to the imperfections on the foil surface as 
well as its curvature and roughness. Due to the non- 
ideality, the data analysis of the experiment assumes that 
the direction of the magnetic moment is an unknown pa- 
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FIG. 4: (color online). Magnetic moment measured by 
MPMS-SQUID for deposition of 3 A Co/ 15 A Pt/3 A Co/ 
30 A Pt on gold foil substrate, including seed layer of 100 A 
Pt deposited at 350C, for in-plane and out-of-plane magneti- 
zation. 



rameter. However the magnitude of the moment can be 
reasonably well controlled as discussed in the following 
subsection. 



D. Magnetic Mass preparation 

The force due to the differing magnetic susceptibility 
of gold and silicon depends on the ambient field Bq which 
magnetizes the drive mass. Bq gets a contribution from 
the Earth's field, but also may be enhanced by the local 
magnetic environment. In particular the magnetic field 
emanating from the ferromagnetic layer on the test mass 
contributes to the drive mass magnetization. The differ- 
ential volume-magnetic-susceptibility of gold and silicon 
Ax is 3 x 10~ 5 , leading to a differential magnetization of 
order AM = j^A\Bq. The correction to the local mag- 
netic field and its gradient near the drive mass surface 
can be estimated by integrating the induced magnetic 
dipole moment per unit volume over the volume of the 
drive mass. Very roughly, the differential magnetic force 
experienced when the magnetic test mass is above gold 
versus above silicon can be approximated as 



AF m c V eS n AM/r 4 



(4) 



where r is the distance to the test mass magnetic film 
from the center of the drive mass bar, m c is the magnetic 
moment on the cantilever, and V c ff is the effective volume 
of the drive mass within range to contribute to the signal. 
To obtain a precise estimate of this force the net direction 
of the magnetic moment on the cantilever must be known; 
otherwise only order-of-magnitude estimates are possible. 
The magnetic measurements described in Refs. @, G3 
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can be explained by Eq. ((4]) to within uncertainties, and 
therefore Eq. ^ can be used to determine the target 
amount of magnetic material on the test mass. 

To obtain a large enough force signal for a reasonable 
drive mass current (10 mA), while also having negligi- 
ble magnetic susceptibility force with the current off, it 
is sufficient to have a magnetic shield attenuation fac- 
tor greater than 100 and a total moment to satisfying 
5 x 10" 15 J/T < to < 5 x 10" 13 J/T. For two layers 
of 3 A Co / 15 A Pt, measurements suggest an area of 
(25 A*m) 2 yields a moment of 5 x 10" 13 J/T. The foil is 
inserted into the FIB and sections of the magnetic layer 
near an edge of the foil are removed with a 100 pA ion 
beam. 12 /im squares of magnetic material are left in 
the locations that will become the centers of the test 
masses and should result in a moment of approximately 
1 x 10" 13 J/T. A high-current 20 nA ion beam is then 
used to cut away the edge of the foil to produce a clean 
edge. The test masses are defined using the high cur- 
rent beam. Care is taken to ensure that the high current 
beam does not directly etch any magnetic material, as 
it has a tendency to deposit material as the beam cuts. 
A thin film of gold covers the magnetic patches on the 
test mass due to this effect. The masses are completely 
free except for a small "tail" . At this point the masses 
are reinserted into the SQUID and re-magnetized out- 
of-plane. This re-magnetization is performed to elimi- 
nate the effects of any demagnetization that may have 
occurred from heating while the high current beam cuts 
through the foil. The resolution of the SQUID is insuf- 
ficient to measure the magnetic moment of a single test 
mass. However, the bulk foil can be seen to attain the 
same value of magnetic moment as originally measured. 
After re-magnetization the masses are mechanically ex- 
tracted from the foil using non-magnetic tungsten probe 
tips. In this way subsequent heating from the ion beam 
can be avoided. The masses are epoxied onto the can- 
tilever using Hysol-Loctite 1C epoxy [16( and again only 
handled with a tungsten probe tip to avoid magnetic con- 
tamination. 



IV. NOISE AND BACKGROUND SOURCES 

Thermal Noise. The minimum detectable force due to 
thermal noise at temperature T is 

\ \k B TBk 
F min = \ , (5) 

where B is the bandwidth of the measurement, k is the 
cantilever spring constant, and /o = loq/(2tt) is its reso- 
nance frequency. This represents a fundamental limit on 
the sensitivity of the measurement technique. 

Vibrational Noise. At the cantilever resonance fre- 
quency, assuming a 4% non-linearity in the bimorph and 
Q of 10 5 , we estimate the required internal vibration 
isolation (between the bimorph and cantilever) to be 



1.3 x 10 5 , employing a simple analytic model [Ic}. Two 
sets of spring-mass stages isolate the platform on which 
the bimorph is mounted from the mount holding the can- 
tilever, providing a factor of 10 6 attenuation at 100 Hz 
and 10 8 attenuation at 300 Hz [7[. The mass stages are 
also connected by thin wires in addition to the springs 
to operate the various sensors and piezo-electric stack. 
If these wires become too taut during the cooling of the 
probe, the vibration isolation system can become com- 
promised. A direct measurement of the vibrational cou- 
pling at 300 Hz is unfortunately not possible due to the 
operating frequency range of the bimorph. The cryostat 
is hung from 1 Hz springs from a thick concrete ceiling 
to attenuate unwanted external vibration at 300 Hz. 

Electrostatics. The stainless steel frame which holds 
the cantilever wafer is grounded to the main probe 
ground, so that together with the metallized shield, a 
Faraday cage surrounds the cantilever. Any voltage dif- 
ference AV built up between the cantilever and shield 
will result in attractive force between them of magnitude 
F = e A(AU) 2 /2ci 2 where A is the area of the cantilever 
surface, and d is the separation between the cantilever 
and shield. Any built up charge on the cantilever can 
produce an electrostatic force. In a measurement using a 
similar shield, the shield motion at 300 Hz was found to 
have an upper bound of approximately 5d ~ 1 pm [H, [TlT | . 
Therefore voltage differences of 100 mV can drive the 
cantilever above thermal noise. Alternatively, if a spuri- 
ous electrical voltage signal Vb s'mojQt with Vb > 100 /.iV 
leaks into the shield, the cantilever can be driven above 
thermal noise. In either case the signal could not be mis- 
taken for a new Yukawa force since it would not vary 
with the spatial periodicity of the drive mass pattern. 

Casimir force. There is a relatively large (~3x 10" 16 
N) Casimir force present at all times between the can- 
tilever and shield membrane. If constant, it does not 
impede a force measurement at the cantilever resonance 
frequency. There is also a Casimir interaction between 
the shield membrane and drive mass ground plane. The 
oscillating bimorph can drive the shield membrane into 
motion, and due to piezo non-linearity, can excite it at 
the cantilever resonance frequency. The shield membrane 
is sufficiently stiff so that the modulation of the Casimir 
force on resonance does not excite the cantilever beyond 
the level of thermal noise ~ 10" 18 N for our measure- 
ment bandwidth [f| @]. Even if the shield were driven 
enough to excite the cantilever on resonance above ther- 
mal noise, the signal would not exhibit the y-equilibrium 
position dependence that a Yukawa signal would have, 
and therefore could not be erroneously interpreted as a 
new force. 

One can also ask whether the shield membrane and 
drive mass ground plane are sufficient to screen the 
direct Casimir force between the drive mass and test 
mass, given the finite conductivity, plasma frequency, and 
thickness of the gold coating. In fact, the drive mass 
ground plane alone is enough to prevent the differential 
Casimir force from the gold and silicon drive mass sec- 
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tions from b eing transmitted at a detectable level. Fol- 
lowing Ref. we employ a reflection-based model for 
computing the Casimir force between two metallic walls. 
We obtain a reduction factor rjp, (denned precisely in 
Ref. which describes the reduction of the Casimir 

force at small separations and for finite conductivity. The 
parameter r/p is defined as Fc = r/pFp where Fp is the 
perfect conductor result. We take a gold test mass of 
thickness 30 fiui as one mirror, and consider the differ- 
ential Casimir force for the situations of the other mirror 
being 25 fiia away and composed either of 0.1 /xm thick- 
ness Au (corresponding to the drive mass ground plane) 
or 100.1 /im thickness Au (corresponding to the drive 
mass ground plane plus gold bar). We find that the dif- 
ference in rjp is 0.0028, corresponding to a differential 
Casimir force of 2.4 x 10 -20 N, allowing a = 1 to be 
probed at A = 20 /im, which is well below the sensitiv- 
ity of the experiment. The shield membrane attenuates 
this Casimir interaction even further, rendering the effect 
negligible. 

Magnetic background. Apart from those discussed in 
Sec. IIII1 there are other mechanisms for a magnetic cou- 
pling that occur with the drive mass spatial periodicity. 
These include the susceptibility-induced interaction be- 
tween the drive mass and the gold part of the test mass 
and eddy currents produced in the drive mass as it oscil- 
lates in the remnant ambient field. Also, any charge built 
up in the silicon sections of the drive mass can generate 
a weak magnetic field as the drive mass moves. All of 
these mechanisms are too weak to produce a measurable 
signal in this experiment. 

Interferometer noise. Sources of electronic and optical 
noise in the interferometer are discussed in Refs. [a. Illj|. 
Electronic noise consists of shot noise from the laser, and 
to a lesser extent Johnson noise from the 10 Mfl feedback 
resistor and fluctuations in the laser current supply, and 
combines to produce typically a few /zV/Hz 1 / 2 . Optical 
noise can be produced by stray reflections in the interfer- 
ometer as well as mechanical shaking of the fiber. Modu- 
lation of the laser at frequencies over 100 MHz serves to 
reduce stray reflections by reducing the coherence length 
of the laser [l8|, [ijj . Mechanical shaking of the fiber can 
mimic cantilever displacement, along with any electrical 
coupling at harmonics of the bimorph drive signal. The 
amplitude of these noise sources is assessed by measur- 
ing with the drive frequency off-resonance from a sub- 
harmonic. However, this voltage noise is not expected 
to vary periodically as the y-equilibrium position of the 
drive mass oscillation is scanned and thus does not mimic 
a Yukawa-like force. 



V. DATA 

The general cooling-down procedure and data acqui- 
sition technique are discussed at length in Refs. [3, Q • 
The resonance frequency fa of the cantilever ranged from 
324.082 to 324.136 Hz while near base temperature dur- 



ing the course of the experimental run which is the cen- 
tral subject of this paper. The base temperature of the 
probe ranged from approximately 11 — 13 K. The quality 
factor as measured using the ring-down technique (see 
for example Ref. [8|) was approximately 85,000. The 
spring constant k of 0.0062 N/m was inferred from the 
measured resonance frequency and calculated mass of the 
test mass. The effective temperature of the cantilever is 
deduced from k and the square of the amplitude spectrum 
near the cantilever thermal peak through the equiparti- 
tion theorem. This effective noise temperature was typi- 
cally measured to be 20 — 25 K at low laser power, about 
10 K above the base probe temperature, as was typically 
seen in previous work Q. The quality factor and effec- 
tive temperature of the cantilever were both found to 
depend on the incident laser power. The laser power was 
reduced below 1 /jW before temperature and quality fac- 
tor stopped changing. The low-temperature fringe height 
(as defined in Table[TJ was correspondingly reduced to ap- 
proximately 250 mV. The voltage noise on the interfer- 
ometer, as assessed by tuning the piezo driving frequency 
100 mHz off resonance, was at the level of 10 -18 N for the 
working fringe height and quality factor. For the gravity 
force measurement, the quality factor of the cantilever 
was adjusted with feedback Q to 9500 — 10000, and the 
effective temperature with feedback was between 2 — 3 
K. The thermal noise limited force sensitivity was ~ 200 
aN/VHi. 

Prior to collecting data in search of non-Newtonian 
gravity (referred to hereafter as 'gravity data'), the ap- 
proximate x- and y-position of the drive mass with re- 
spect to the test mass is coarsely determined by using the 
magnetic force signal, and it is verified that the test mass 
is positioned over the central meandering region of the 
drive mass. Gravity data are then collected as a function 
of the y-equilibrium position of the drive mass oscillation 
for a given z-separation. For each point of gravity data 
collected, a magnetic force measurement is done at the 
same location by applying current across the drive mass. 
Then the magnetic force is scanned in the y-direction to 
determine the y-position of the closest magnetic mini- 
mum. In this way the relative y-position between grav- 
ity data points can be determined even in the presence 
of drifts in the y-capacitance sensor reading of a given 
y-position that can occur over the time scale of days, as 
demonstrated in earlier experimental runs with a nickel 
magnetic layer on the test mass Q. Having the mag- 
netic calibration (switchable in situ) thus permits mul- 
tiple days of data to be analyzed together, a significant 
improvement over previous renditions of the experiment. 

The data set is shown in Fig. [5l The magnitude and 
phase of the magnetic (drive mass current on) and pos- 
sible Yukawa (drive mass current off) forces at the third 
harmonic are studied as the y-equilibrium position of the 
drive mass oscillation is varied. The 17 data records 
each consist of 24 — 120 minutes of averaged data. The 
amount of averaging time chosen varies considerably due 
to the limited duty cycle in the experiment. Environ- 
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y-equilibrium position (|im) 



FIG. 5: (color online). Main data set showing statistical and 
systematic uncertainties in the measured force. Diamonds 
show signal with current off. Circles show signal with 10 mA 
current in the drive mass. Lines are guides to the eye. Er- 
ror determination is described in Sec. IVII The ^-coordinates 
shown for the data set on the horizontal axis are the displace- 
ments relative to one particular magnetic force minimum cho- 
sen as y = 0. We note that these coordinates vary by an un- 
known offset with respect to those shown in the calculations 
in Figs. [2]and[l 



mental noise and possibly the internal (bimorph-related) 
vibrational noise varied considerably over the course of 
the data run. By 'environmental noise' we here refer to 
noise that is not associated with the motion of the bi- 
morph actuator. In particular, environmental noise may 
have come from the several nearby large scale construc- 
tion projects that were underway during various stages of 
the data run. As a test of the environmental noise, mea- 
surements were taken with the bimorph-actuator turned 
off. Internal vibrational noise can in principle be caused 
by a short of the vibration isolation springs due to taut 
probe wiring. Although considerations are taken at room 
temperature when the system is open to prevent taut 
wiring, it can be difficult to ensure that no vibration 
isolation shorts develop after cooling. The amount of 
internal (bimorph-related) vibration noise can therefore 
vary based on the x, y, and ^-separation and tilts of the 
bimorph stage with respect to the mass-stage holding 
the cantilever. The effects of vibration due to internal 
isolation shorts were examined with the test and drive 
masses far (~ 1 mm) withdrawn from each other and the 
bimorph moving. These noise sources limited the amount 
of useful averaging time over the experimental run. The 
data in Fig. [5] were taken only at time periods when such 
noise was minimal. Since we are not searching for rare- 
event physics, it is legitimate to discard those other data 
for which large disturbances were clearly present on the 
system. 

The phase of the magnetic force undergoes a change 



of 7r at the magnetic force minimum at y — 0. If a pure 
Yukawa signal were present, one expects to see two phase 
changes of it in the drive mass current-off (abbreviated 
hereafter as 'current-off') data for every phase change 
in the magnetic (current-on) signal, and two force min- 
ima for every magnetic force minima. The y-extent of 
the data set should include three Yukawa phase changes, 
each separated by 100 /mi in the presence of such a force. 
The current-off data therefore shows no clear Yukawa 
signal in Fig. [5] In fact, the phase of the current-off 
data remains relatively constant over the data set. In 
the presence of thermal noise alone, the phase would be 
expected to be distributed randomly. The bound on any 
Yukawa signal associated with the data set is presented 
in Sec. IVII The strength of the magnetic moment can 
be roughly estimated from the current-on response. Sim- 
ulation results indicate a magnetic moment of approxi- 
mately 1.3 — 3.7 x 10~ 14 J/T, depending on the direction 
of the moment. 

The relatively constant phase of the current-off data 
could have a number of possible origins. The driving sig- 
nal applied to the piezo actuator could have some small 
electronic leakage to the piezo near the cantilever. The 
internal vibrational noise due to the bimorph oscillation 
can tend to occur at a common phase given a common 
mechanism. Bimorph non-linearity is also responsible for 
some shield motion at the cantilever resonance frequency. 
Although the associated modulation of the Casimir force 
is too small to be observed, excess charge on the can- 
tilever could couple to small displacements of the shield 
membrane driven by the bimorph. The patch effect can 
produce local potential variations in the gold coating on 
the shield membrane which can vary with respect to the 
potential on the gold test mass [2(|. In any case, the 
expected 100 fim spatial periodicity of a Yukawa-like sig- 
nal can be distinguished from a background force with a 
relatively constant phase. Therefore, an accurate bound 
on a Yukawa-like interaction can be determined even in 
the presence of such a background force, as discussed in 
Sec. EH 



VI. ERROR ANALYSIS 

The experimental error in the determination of cither 
a bound on or discovery of a Yukawa-type correction to 
Newtonian gravity can be sorted into two types. First 
there are the statistical and systematic uncertainties in 
the measurements of the force -Fmoas ■ The second type is 
the systematic error which enters the finite-element cal- 
culation of the expected Yukawa force signal Fy- We 
treat the experimental uncertainties in a Monte-Carlo 
fashion, based on the method described in Ref. [7j. In 
order to obtain the function a(A) for a Yukawa force that 
best fits the measured data we consider a least-squares 
fit of the real and imaginary parts of the data to a set 
of calculated forces in the Monte-Carlo simulation. The 
four fitting parameters, as a function of A, are the re- 
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TABLE I: Uncertainty in measured values of forces. The 
voltage on the interferometer is given by V^ nt f(d) = V c — 
-^jp- cos(47rd/A), where d is the distance between test mass 
and the end of the optical fiber, and A is the laser wavelength. 
The fringe height V pp = (Vm ax — Knin) and the fringe center 
V c is the mean of V max and V m i n . The fringe position is the 
deviation from the fringe center. 



Parameter 



Value 



Error Units dF% 



Fringe height 
Fringe position 
fo 

Q 
k 

Fiber alignment 



.228-.472 


324.082-324.136 
9500-10100 
.0062 




002-.004 
3 
.002 
400 
.0003 
22 



V 
% 
Hz 

N/m 
/im 



1 

0.2 
0.7 
4.1 
4.5 
12 



TABLE II: Uncertainty in tilt and z-separation. Total tilts 
of the drive mass relative to the shield plane are expressed in 
vertical distance over the 1.8 mm length or 1.3 mm width of 
the drive mass. 



Parameter 



Value la Error Units 



z-separation 

Shield-to-drive mass 11-15 1 /im 

Edge-edge between masses 27-31 3 /jm 

Mean edge-edge separation 29 3 /im 

MC Input 29 3 /im 
Relative to the shield plane 

Total tilt in y-z plane 2 5 /im/length 

Total tilt in x-z plane 3 5 /im/width 



Total 



13 



measurement added in quadrature is approximately 13%. 



suit a, an offset variable yo to account for the lack of 
absolute y-position information between the drive mass 
pattern and test mass due to the unknown direction of 
the magnetic moment, and the real and imaginary parts 
of a constant offset force (Ro, Io). It is reasonable to in- 
clude a constant offset force to account for vibrational or 
other-type backgrounds that do not depend periodically 
on the y-equilibrium position of the drive mass. 



A. Uncertainty in the measured force 

Statistical errors. The statistical error on each data 
point comes primarily from thermal noise at the effective 
cantilever temperature and is the most significant contri- 
bution (90%+) to the total la error bars shown in Fig. [5] 
The variation in the standard error reflects the amount of 
averaging time for each data point, which varies from 24 
to 120 minutes. The statistical error decreases with the 
square root of the numbers of samples. There is a natu- 
ral time scale for the duration of statistically dependent 
samples given by the ring-down time of the cantilever 
t = nf Q. We divide the data records into subsections 
of length 30 seconds ~ 3r and study the cross correla- 
tion between these subsections to verify their statistical 
independence, by the method discussed in Ref. 0]. 

Systematic errors in the measured force. The system- 
atic errors in the measured force are listed in Table U 
The error is dominated by the uncertainty in the exact 
position of the fiber over the top surface of the test mass. 
Considering the cantilever mode shape we estimate the 
corresponding la error in the measured force as 12%. 
The uncertainty including other systematics in the force 



B. Uncertainty in the calculated force 

Uncertainty in vertical separation and tilt. The cal- 
culated Yukawa force strongly depends on the vertical 
separation between the drive and test masses due to its 
exponential form. The tilt of the drive mass with re- 
spect to the shield wafer will also effectively increase the 
z-separation, by an amount that depends on the x and y 
position of the test mass with respect to the drive mass 
edges. Table |n] indicates the experimental uncertainty in 
the z-separation and tilts. The tilts and z-separation are 
determined using methods described in previous work [tJ . 
From this table, inputs to the Monte-Carlo are derived. 
The Monte-Carlo input for z-separation is the mean edge- 
edge separation, adjusted to account for the tilt in the 
drive mass silicon plane jlQ |. Variations in z-separation 
between data points are very similar to variations within 
a particular data point (1 /im) and typically are less than 
the overall uncertainty in z-separation (3 /im). Also for 
the tilts, variation of tilts across data points are much less 
than the overall uncertainty in the tilt. Therefore for the 
Monte-Carlo we can treat the z-separation and tilts of 
the entire data set as their mean value and add a Gaus- 
sian random error with their standard deviation. This is 
a reasonable approximation since the standard deviation 
in the measured tilt or z-separation over the 17 point 
data set is significantly less than the uncertainty in the 
determination of the tilt or z-separation at any one data 
point. A more detailed accounting of the determination 
of tilt and separation uncertainty appears in Ref. [13]. 

Uncertainty in bimorph amplitude. The mean bimorph 
amplitude recorded for the data set was 119 /tm with a 
standard deviation of 1.3 /im. The la uncertainty due 
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TABLE III: Errors due to masses 



Parameter 



da% 



Volume of Test Mass 4.5 
Voids in Drive Mass 2.5 
Drive mass polished Au/Si boundary 1 



Uncertainty in the expected phase of a Yukawa-type sig- 
nal. Simulation indicates that the phase of any Yukawa 
signal will be identical to the phase of the magnetic force 
modulo 7r. The uncertainty in the expected phase (mod- 
ulo 7r) of a Yukawa signal is determined from the uncer- 
tainty (0.1 rad) in the measured phase of the magnetic 
calibration signal. 



C. Summary of inputs to the Monte-Carlo 



TABLE IV: Uncertainty in Masses and Geometry 



Parameter 


Value 


Error 


Units 


Test Mass length 


53 


1 


fxm 


Test Mass width 


54 


1 


fim 


Test Mass height 


27 


1 


/im 


Volume of Test Mass 


77000 


3500 


/im 3 


Density of gold 


19300 





kg/m 3 


Total mass of test mass 


1.50 


0.07 





to the motion calibration is 10 /im. The method for cal- 
ibrating the bimorph and for determining the error in 
bimorph amplitude is discussed at length in Ref. Q. 

Uncertainty in geometry and density of masses. The 
uncertainty in the geometry of the drive mass is discussed 
at length in Ref. |8j. The uncertainty in the test mass 
geometry is determined by a combination of optical and 
SEM imagery. Table [TTT1 lists their relative contribution 
to the uncertainty in a, which is included as a multiplica- 
tive scaling error to the best fit results, using the method 
discussed in Ref. 0- The effective length, width, and 
height shown in Table IIVI are determined by considering 
the approximately 84° sloping profile of the test mass 
walls measured by the SEM and accounting for the ex- 
tra mass present near the final-release corner. A SEM 
micrograph showing the mass attached to the cantilever 
after being used in the data set is shown in the inset of 
Fig.ffl 

Uncertainty in relative y-position of data points. For 
the Monte-Carlo analysis to be optimally constraining, 
the relative y-separation of the data points must be 
known to within a few microns. The y-coordinate of each 
data point has a la uncertainty of 2.5 /im. This accounts 
for the uncertainty in the y-position of the magnetic min- 
imum, as well as the uncertainty in the y-position at the 
data collection point. 



The set of geometrical inputs (z-separation and tilts in 
the y-z and x-z planes, respectively) for the Monte-Carlo 
calculation are summarized in Table [TXJ For each of the 
values of A considered (4,6,10,18,34,66 ^m), the Yukawa 
force versus y-displacement is calculated for 320 elements 
of a Gaussian distribution of these geometrical parame- 
ters, with the means and standard deviations as indi- 
cated. For the mean value of bimorph amplitude, the cor- 
responding Zuj force versus y-displaccmcnt is determined. 
For each of these curves, a least-squares fit to the mea- 
sured forces and y-positions is performed, yielding the 
best-fit a, yo, Ro, and Iq. The measured force, measured 
y-position of the data points, and expected phase of the 
Yukawa type signal arc then 'dithered' (with a Gaussian 
random value added) 200 times according to the Gaus- 
sian distribution of their statistical and systematic errors 
as described in the previous subsections. For the uncer- 
tainty in phase, an element from a Gaussian distribution 
is randomly assigned to the expected phase (modulo it) of 
a Yukawa type signal, according to the uncertainty in the 
phase of the measured magnetic signal. Each dithered set 
of data is then fit to yield a new best-fit value of a. This 
results in a set of 320 x 201 = 64,320 least-squares fits 
of the data to the calculated curves, producing 64,320 
best- fit values of a for each A. Finally, these values of a 
are appropriately scaled to account for the uncertainty 
in the bimorph amplitude and the uncertainty in the vol- 
ume of the drive and test masses, by using 50 samples 
taken from Gaussian distributions representing these un- 
certainties. The treatment of these multiplicative errors 
is further discussed in Ref. [l(|. The end result is a his- 
togram of 3,216,000 best fit a for each of the six values 
of A explicitly considered. The numbers of samples are 
chosen to be sufficiently large so that the means and 95- 
th percentiles of the histograms do not vary significantly 
with the numbers of samples. 



VII. RESULTS 

The histogram resulting from the Monte-Carlo analysis 
for A = 10 /im is shown as an example in Fig. [6] To illus- 
trate a typical fit of the force calculation to the data we 
consider an element near the mean of the histogram for 
A of 10 [im, corresponding to a — 4450. Figure [7] shows 
the real and imaginary parts of the data along with sta- 
tistical and systematic error bars. Also shown is the best 
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FIG. 6: (color online). Histogram of best-fit a results for 



A — 10 /im. 
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FIG. 7: (color online). Real and imaginary parts of measured 
force. A typical Monte-Carlo fit result (best fit Yukawa signal 
plus offset) for A = 10 ^m, a = 4450, is shown as a solid line. 
Error bars (la) on the measured forces are also shown. 



fit Monte-Carlo result for comparison. Figure [S] shows 
the amplitude and phase of the data along with statisti- 
cal and systematic error bars. The solid line shows the 
best fit Monte-Carlo result to the data. We note that the 
actual fit is performed on the real and imaginary parts of 
the data — the amplitude and phase are plotted for illus- 
tration. For reference, the dotted line indicates the best 
fit result after removing the constant offset force. The 
fit is performed with the offset included and the offset 
is subtracted afterwards. The data set is clearly best fit 
with a constant offset force included, as evidenced by the 
concentration of phase near 1 rad for the entire data set. 
The fact that the best-fit phase does not vary much re- 
flects the adding of the best-fit Yukawa component (with 
an amplitude of 0.5 aN) to the larger background best-fit 
offset rms force vector (2.1 aN, 3.1 aN). 

The mean values of a and 95% confidence exclusion 
bounds for Yukawa type forces are derived from the set 
of six histograms and listed in Table [V] The ratio of 
the 95% confidence limit to the mean value increases 
as A decreases, reflecting the increased significance of 
the uncertainties in z-separation and tilt. The means 
of each histogram represent the most likely value of any 
Yukawa force that may be present. Upcoming experi- 
ments will exhaustively re-examine this parameter space. 
The results represent an improvement of the constraints 
on Yukawa forces by approximately half an order-of- 
magnitude over our previously published work Q- We 
also note that the mean values of a from the histograms 
in Ref. Q are now excluded at better than the 95% level 
by the present work. Our final results are shown in Fig. 
O We show only the 95% confidence exclusion results 
from Table El 
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FIG. 8: (color online). Amplitude and phase of measured 
force. A typical Monte-Carlo fit result (best fit Yukawa signal 
plus offset) for A = 10 fim, a = 4450, is shown as a solid line. 
Error bars (ler) on the measured forces are also shown. Dotted 
line shows fit result with best-fit offset-force subtracted. The 
amplitude of the best-fit offset-force is significantly larger than 
the amplitude of the Yukawa component. 



VIII. FUTURE DIRECTIONS 

One limiting factor of the experimental run had to do 
with the duty cycle of the experiment. Improvements to 
the duty cycle could be made by adding more automa- 
tion to the experiment or transferring the apparatus to 
a generally less noisy environment. A substantial im- 
provement could involve the redesign of the cantilever 
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TABLE V: Experimental limits on Yukawa forces 



X(fim) Mean(MC) a 95% exclusion a 



4 


8.6 x 10 6 


3.1 x 10 7 


6 


1.6 x 10 5 


4.6 x 10 s 


10 


5.6 x 10 3 


1.4 x 10 4 


IS 


5.1 x 10 2 


1.1 x 10 3 


34 


1.2 x 10 2 


2.5 x 10 2 


66 


7.0 x 10 1 


1.5 x 10 2 




X (|jm) 



FIG. 9: (color online). New limits on Yukawa-forces. The 
shaded area refered to as 'Stanford 3' represents the 95% ex- 
clusion region yielded with the present data set, listed in Ta- 
ble |V] Experimental limits derived in this work are on \a\. 
The previous results labeled 'Stanford 1' and 'Stanford 2' are 
presented in references Hj] and 0], respectively. Also shown 
are theoretical predictions and other exp erimental results as 
indicated in Refs. [U, [H [H, HI El Hf . 



to allow a larger area test mass. Such an improvement 
is underway in a second-generation rotary-drive experi- 
ment [27j , expected to be one to two orders of magnitude 
more sensitive than the present apparatus. A number of 



other improvements could be implemented resulting in 
more marginal gains. A doubly-clamped cantilever beam 
loaded with a test mass in its center could reduce the sys- 
tematic error associated with the transverse position of 
the fiber above the test mass. In this configuration, any 
transverse displacement of the fiber along the test mass 
length affects the measured vertical displacement to a 
lesser degree, as the cantilever mode shape is quadratic 
rather than linear at the mass-attachment point. An 
additional piezo-electric transducer that, unlike the bi- 
morph, can be operated at the cantilever resonance fre- 
quency could be installed near the bimorph to further 
characterize the effectiveness of the vibration isolation. 
Voltage noise on the interferometer could potentially be 
reduced by improving the vibration isolation between the 
optical fiber in the cryostat and the bimorph actuator. 
A new set of shield wafers could be fabricated from a 
double-polished silicon wafer to allow a sharper reflec- 
tion of the drive mass image during the room tempera- 
ture alignment procedure. This could potentially reduce 
uncertainty in the tilt, resulting in decreased effective 
separation of the masses. 



IX. CONCLUSION 

To date, the gravitational interaction remains as one of 
the least well-understood and least-tested aspects of fun- 
damental physics. In this paper we have presented the 
latest results from the first-generation Stanford micro- 
cantilever experiment, which in sum has improved the 
limits on new Yukawa-type forces at 20 /im by over four 
orders of magnitude since its conception. Our most re- 
cent experimental constraints on Yukawa-type deviations 
from Newtonian gravity represent the best bound in the 
range of 5 — 15 /j,m, with a 95% confidence exclusion of 
forces with \a\ > 14,000 at A = 10 ^m. 
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